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COMPLEX EXCEPTIONAL ORTHOGONAL POLYNOMIALS 
AND QUASI-INVARIANCE 

WILLIAM A. HAESE-HILL, MARTIN A. HALLNAS, AND ALEXANDER P. VESELOV 


Abstract. Consider the Wronskians of the classical Hermite polynomials 
■.= Wr{H,{x),Hk,{x)...,Hk^{x)), Z S Z>o, 

where ki = Aj + n — i = 1,... ,n and A = (Ai,..., An) is a partition. 
Gomez-Ullate et al showed that for a special class of partitions the correspond¬ 
ing polynomials are orthogonal and dense among all polynomials with certain 
inner product, but in contrast to the usual case have some degrees missing 
(so called exceptional orthogonal polynomials). We generalise their results to 
all partitions by considering complex contours of integration and non-positive 
Hermitian products. The corresponding polynomials are orthogonal and dense 
in a finite-codimensional subspace of C[x] satisfying certain quasi-invariance 
conditions. A Laurent version of exceptional orthogonal polynomials, related 
to monodromy-free trigonometric Schrodinger operators, is also presented. 


1. Introduction 


Consider polynomials Pn{x) G ]R[a:] of degrees n = 0, > satisfying the orthog¬ 

onality relation 

{Pm^ Pn') — ^mnQnj 

where the inner product of polynomials is defined by a real integral 

(1) {p,q)-= p{x)q{x)w{x)dx 

J a 


for some positive weight function w. Suppose that there exists a second order dif¬ 
ferential operator 

having these polynomials as eigenvectors: 


Tpn{x) = EnPn{x), n = 0, 1, • • ■ • 

A classical result due to Bochner [3] says that in that case the sequence of polyno¬ 
mials Pnix), n G Z>o, must coincide (up to a linear change of x) with one of the 
systems of classical orthogonal polynomials of Hermite, Laguerre or Jacobi. 

Gomez-Ullate, Kamran and Milson m considered the following variation of 
Bochner’s question. Let us assume now that in the previous considerations n 
belongs to a certain proper subset S C Z>o such that Z>o \ S' is finite. To 
make this non-trivial they added the following density condition: the linear span 
U = (pn : n G S) of the corresponding polynomials must be dense in K[a;] in 
the sense that if (p,p„) = 0 for all n G S then p = 0. In that case the sequence 
Pnix), n G S is called a system of exceptional orthogonal polynomials. 
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The main example of such polynomials are exceptional Hermite polynomials [12] 
having the Wronskian form 

(2) HxAx)--=WTiHi{x),Hk,ix)...,HkAx)), ^ € Z>o, 

where Hi{x) are classical Hermite polynomials, A = (Ai,..., A„) is a double parti¬ 
tion and 

ki = Xi + n — i, f = 1,..., n. 

The double partitions have the very special form 

where fj, = (/ri, /r 2 ,..., p-fe) is another partition with n = 2k (see [S]). According to 
Krein and Adler [T] this guarantees that the corresponding Wronskian 

(3) Wx{x)=WviHk,ix)...,HkJx)) 

has no zeroes on the real line and thus determines a non-singular weight function 

(4) w{x) = W^^{x)e~^ . 

The geometry of the complex zeroes of the corresponding Wronskians is quite in¬ 
teresting and was studied by Felder et al. in [8]. 

One of the goals of our paper is to find a proper interpretation of the exceptional 
Hermite polynomials (§ for all partitions A. As we will see, this will naturally lead 
us to the notion of quasi-invariance, which appeared in the theory of monodromy- 
free Schrodinger operators, going back to Picard and Darboux and more recently 
revisited by Duistermaat and Griinbaum [7]. In certain classes such operators 
were explicitly described in terms of Wronskians in [3 H] [Ml [10] . Grinevich and 
Novikov studied the spectral properties of these and more general singular finite- 
gap operators and emphasized the important link with the theory of Pontrjagin 
spaces (see m and references therein). Our paper can be considered as dealing 
with the implications of all these results for the theory of exceptional orthogonal 
polynomials. 

More precisely, we first complexify the picture by considering the vector space 
V = C[z] and replace the inner product Q by a Hermitian product of the form 

(5) {p,q)--= / p{z)q{z)w{z)dz, 

Jc 

where q{z) := q{z) is the Schwarz conjugate of the polynomial q{z), C C C is 
a contour in the complex domain and w{z) is a complex weight function. The 
condition that this product is Hermitian implies certain restrictions on the contour 
C and function w(z) (see section 2). It also requires certain restrictions on the 
set of polynomials for which the product is well defined. As it turned out, such 
polynomials form a subspace U C V of finite codimension defined by some quasi¬ 
invariance conditions. Similarly to m we say that the polynomials Pniz:), n G S, 
form a system of complex exceptional orthogonal polynomials if their linear span is a 
subspace of U that is dense in U in the sense that {p,Pn) = 0 for all n G S, implies 
that p = 0. 

We will show that the Wronskians ([^ satisfy this criteria for every partition A 
and a suitable choice of C with w given by Q . For a double partition A we can take 
as a contour C the real line with U = V and recover the results by Gomez-Ullate 
et al [T^ . 
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Note that the corresponding Hermitian form is positive definite only for double 
partitions, otherwise we always have polynomials with negative norms. The appear¬ 
ance of negative norms for singular potentials was first emphasized by Grinevich 
and Novikov m- 

We also consider the Laurent version of our approach. Some Laurent versions of 
orthogonal polynomials are already known in the literature (see e.g. [6] and refer¬ 
ences therein), but our approach is different since it is not based on Gram-Schmidt 
procedure. Similarly, it does not fit into the theory of orthogonal polynomials on 
the unit circle initiated by Szego who considered the case of usual polynomials. 

Gonsider the Laurent polynomials A = C.[z,z~^] and the following complex bi¬ 
linear form on A: 

where G = {z S C : [z] = 1} is the unit circle. The standard basis z", n S Z, 
satisfies the Laurent orthogonality relation 

{z’",z^) = dk+i,o, k,lGZ. 

We consider more general forms 

( 6 ) 

where is the circle defined by \z\ = fi and w{z) = W{z)~‘^, with W{z) some 
Laurent polynomial. For this form to be well-defined, we need to assume that P, Q 
belong to a suitable subspace of quasi-invariants ^ C A of finite codimension. 

Let a;; be a finite subset of N. Suppose that P„ G A, n G Z, satisfy the Laurent 
orthogonality relation 


(7) 


{,Pk ) Pi ) — ) k^ I ^ Zj 


but Pn is proportional to P-n for n € 1C, which implies that the corresponding 
hn = 0, and thus Pn is orthogonal to all Pk, k G Z. If the minimal complex 
Euclidean extension of the linear span of P„, n G Z, coincides with the subspace of 
quasi-invariants then we call them exeeptional Laurent orthogonal polynomials. 
The need to consider such an extension is the novelty of the Laurent case, which is 
related to the fact that the corresponding form is degenerate on the linear span of 
Pnt n G Z. 

We present an example of such polynomials corresponding to the trigonometric 
monodromy-free Schrodinger operators |3]. Namely, for any set k = {fci,...,A:„} 
of distinct natural numbers ki > k 2 > ■ ■ ■ > kn > 0 and any choice of complex 
parameters a = (oi,..., a^), afc G C \ {0}, we define the Laurent polynomials 


( 8 ) P.,a-A^) = 


^’/ci(ai;z) 


D^k2{a2;z) 


D^k 


(a„;z) z' 
^(a„;z) Pz' 


P’^$fc,(ai;z) P’^$fc,(a2;z) ••• P-$fc„(a„;z) D^z^ 

where $fc(a; z) = az’^ + a~^z~^, fc G N and D = z-^. 

When parameters ak satisfy the condition |afe| = 1 for all A: = 1,... n we intro¬ 
duce a Hermitian form on a certain subspace of quasi-invariant Laurent polynomials 
^K,c Slid show that the minimal Hermitian extension of the linear span of P^.a-.h 
I G Z, coincides with the subspace of quasi-invariants and is dense in 




4 


W. A. HAESE-HILL, M. A. HALLNAS, AND A. P. VESELOV 


2. Complex exceptional Hermite polynomials 


In this section we consider the polynomials as defined in for general 

partitions, i.e., we do not require that A is a double partition. We shall refer to 
these polynomials as complex exceptional Hermite polynomial or CEHPs for short. 
Although these polynomials have real coefficients for the partitions which are not 
double it is natural to consider them as elements of the complex Hermitian vector 
space because the contour of integration in (§ is complex. 

We begin by recalling how H\y are obtained by a sequence of Darboux trans¬ 
formations from the classical Hermite polynomial Hi. The starting point is the 
classical fact that the functions 

(9) ^i{z)=Hi{z)e-^"/\ lGZ>o 

have the eigenfunction property 

+ z'^tp = {21+ z G C, 

and satisfy the boundary conditions 


lim 

Re 2—>-±oo 


tpi{z) = 0 


(for any fixed value of Imz). We will choose the normalisation of Hermite poly¬ 
nomials such that the highest coefficient of Hi{z) is 2^ and all the coefficients are 
integer: 

Ho = 1 , Hi= 2z, H 2 = 4z2 -2, H 3 = 8z^ - 12z, Hi = 16z'‘ - A8z^ + 12 ,.... 


As is well known, after n consecutive Darboux transformations at the levels 
kn < kn-i < ■ ■ ■ < ki, where ki = \i + n — i, i = 1,... ,n, one arrives at the 
Schrodinger operator 

(10) - 2^(logWr('(/'fei,.. +^^ 

which satisfies the intertwining relation 

o ^ = =S^A o ^A, 


where the intertwining operator acts according to 

Wr(V',V'fei,---,V'/c„) 


( 11 ) 


^aV' = 


Wr(i/>fci,...,-0fe„) 

see e.g. mm- It follows that the functions 


( 12 ) 




Wr(Vifcj,.. . ,'i/'fe„) 

have the eigenfunction property 

^aiAa.z = (2Z -I- l)'0A.i- 

Substituting (|^ in (12) and using the general property Wr((//i,..., <;/„) = 
(;"Wr(/i,..., /„), one finds that 

— x^j2 

(13) i/’A,; = H\,i , 

where Hx i are given by §• 
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O X o - •-O • O ► E 

I 3 5 7 9 II 13 


Figure 1. The eigenvalues of the first few eigenfunctions for A = (1). 


By a direct computation, it is readily inferred that H\ i is an eigenfunction of 
the operator 


T =-^ 
^ dz^ 




WxJ dz 


Wx „ W( 
^ +2z—^ 


Wx 


Wx 


with eigenvalue 2/ + 1. 


Example 1. Consider the special case A = (1), which corresponds to the Schrodinger 
operator 




.—’/a 


dz^ 

d?_ 

dz"^ 


+ z^ + 


)) 


+ 2. 


Already in this simple example, we obtain eigenfunctions (121 with a singularity 


on the real line (at z = 0^. Indeed, this can be seen explicitly by writing out the 
first exceptional Hermite polynomials H{i)^k = Wr('0fc,'(/ ji) and the corresponding 
few eigenfunctions : 




V'( 1).0 



^(i),2 — ~(2 + 4^^), 

= -16^3, 


Fr(i),4 = 12(1 +4^2 _4 z4), 
64^3(5-2^2), 

g = 40(3 + 18 z 2 - SGz'i + Sz®), 


V’(l).6 = 


V'(i),2 = -^e-V2, 

'0(1).3 = -Ifiz^e-^'/z, 

V,(4),4 = H(l±4^e-.V2, 
^(1).5 = 64^2(5 - 2z2)e-^V2, 

_ 40(3-H8z^-36z"^-L8z°) -z^/2 


More generally, using the fact that 

WT{'lp2l+l,1pl)i-z) = -WT{tlj2l + l,i’l)iz), I e N, 

as well as the fact that each classical Hermite polynomial H 2 i{z), I € Z>o, has a 
nonzero constant term, it is readily seen that '0(i),;(a;) is regular on the whole real 
line if and only if I S Z>o\{l} is odd. The eigenvalues of the first few eigenfunctions 
are given in Figure^ where open and filled circles indicate that the corresponding 
eigenfunctions are singular and non-singular, respectively. In addition, the cross 
represents the eigenvalue removed by the Darboux transformation. 

Note that in the theory of quantum Calogero-Moser systems (of which this exam¬ 
ple is the simplest case) only non-singular solutions are considered (see e.g. m)- 
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We will now use the fact that is obtained as the composition of first order 
intertwining operators. To be more specific, let us introduce the short hand notation 

Wm = Wr('(/)fc^, .. = Wr('(/','0fe„, • ■ ■,'4’kJ, 

(where it is convenient to allow m = n + 1 and set Wn+i = 1, Wn+iitp) = tp), and 
recall the standard identity 

“ l^m-1 (V') TO > 1. 

ax ax 

Then it is readily verified that 

( 14 ) = DiO ■ ■ ■ O DmO ■ ■ ■ o Dn, 


with 

(15) 




y _ y / -^ra \ 

dz dz \ Wra+l ) 


For our purposes, a key notion is that of trivial monodromy, see e.g. [18) . A 
Schrodinger operator = —d^/dz"^ + u(z), whose potential u is a meromorphic 
function of z, is said to have trivial monodromy if all solutions of its eigenvalue 
equation 


(16) 


JP’pj{z) = E'ip{z) 


are meromorphic in z for all E. 

We recall that every monodromy-free Schrodinger operator ^ with a quadrati- 
cally increasing rational potential is of the form (10) for some partition A. The fact 
that each Schrodinger operator has trivial monodromy is easily seen. Indeed, 
in the special case u(z) = z^ all eigenfunctions are entire, and trivial monodromy 
is preserved under (rational) Darboux transformations. The converse result is due 
to Oblomkov M- 

Duistermaat and Griinbaum [7] obtained local conditions for trivial monodromy. 
Specifically, in a neighbourhood of a pole z = Zi the potential u(z) must have a 
Laurent series expansion of the form 


4 ^) = Cr{z - Z,Y 


r>-2 


with 

and 


c _2 = mi(mi + 1) for some to^ € N, 


C2j_l=0, Vj = 0, 1, . . . ,TOi. 

In addition, every eigenfunction tp has a Laurent series expansion of the form 

OO 

yz) = (z-zO-™“E^-(^-^*r, 

r=0 


with 

d2j-l=0, Vj = I,...,TOj. 

We proceed to consider the implications for the CEHPs H\ i. Let Z\ be the set 
of zeros Zi G C of the Wronskian Wx{z) with multiplicities mi G N. In addition, 
we need the subset Zf C Zx obtained by restriction to Zi G K. We say that a 
meromorphic function is quasi-invariant at the point z = Zi with multiplicity 
TOi if it satisfies the following two conditions: 
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(1) tp{z){z — ZiY^' is analytic at z = Zi, 

(2) {Y{z){z - for all j = 

The second condition can be rewritten as 

YiaYz)) = (-ir'Y(z) + Oiiz - z^% 

where (Ji(z) = 2zi — z is the reflection with respect to Zi. This explains the termi¬ 
nology. 

Introducing the subspace 

Qa = < P G C[2;] : Y(z) := is 

it follows from the above that the C-linear span 

U\ = {Hx i : I G Z>o \ {fci,..., kn}} 

belongs to Qx- From Proposition 5.3 in [T^], we recall that the codimension of Ux 
in C[z] is equal to |A|. On the other hand, |A| is the degree of Wx{z), and therefore 
the number of quasi-invariance conditions that any p G Qx should satisfy. This 
yields the converse inclusion, and thus the following result. 

Proposition 1. The C-linear span of CEHPs coincides with polynomial quasi¬ 
invariants: 

Ux = Qx- 


quasi-invariant at z = Zifizi G Zx 


Whenever A is not a double partition, the Wronskian Wx{z) will have one or 
more real zeros [T] , so that the weight function is no longer non-singular on the 
real line. To resolve this problem, we replace the standard contour K by a shifted 
contour C = *^ + M and consider a corresponding Hermitian product ([^. As will 
become clear below, to ensure that the product is Hermitian we need to restrict 
attention to the following subspace of quasi-invariant polynomials: 


Qa,r = {p & C[z] : ip{z) := p{z)- 




is quasi-invariant at z = ZifiZi G Z^ > . 


Wx(z) 

By counting quasi-invariance conditions, we obtain the next proposition. 


Proposition 2. The codimension of Qx in Qa,r is |A| — Ylz 
We are now ready for the main definition of this section. 


Definition 1. Let f gM. be such that 

(17) 0 < lei < |Im z,|, \Jz,GZx\Zl. 

Then, we define a sesquilinear product (•,•) on Qa,k setting 

f 

(18) {P,q)= Piz:)qiz)---^^dz, p,qGQx,R- 

Now we will show that the product does not depend on the specific choice of f. 
We find it worth stressing that this important property relies on our restriction to 
the subspace Qa,k- 

Proposition 3. For any p,q G Qa.r, the value of {p,q) is independent of f G M. 
provided the condition 0 is satisfied. 
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Proof. Let denote the integral in the right-hand side of (181. By Cauchy’s 
theorem, it suffices to show that = 0 for some ^ satisfying 0- From 

the residue theorem, we deduce that the difference between the two integrals is 
proportional to 

\ AV 

We claim that each of these residues vanish. In fact, we have the following more 
general result. 

Lemma 1. If ip, (p are quasi-invariant at z = Zi with multiplicity rrii, then 

Res (jp{z)(p{z)) = 0. 

Indeed, it follows from Condition (2) above that 

{ip{z)(p{z){z - I 


3=0 


= 0 . 


□ 


It is now straightforward to show that Definition yields a Hermitian product. 
Proposition 4. The sesquilinear product (•, •) is Hermitian: 


{p,q) = {q,p), Vp,9gQa.r- 

Proof. In what follows, we find it convenient to use the notation 


,{z) = 


e-P/2 

wfily 


and use a subscript to indicate the choice of ^ in ( |18[ ). Since the classical Hermite 
polynomials have real coefficients, it is evident from (§ that w{z) = w{z). Hence, 
we have the following equalities: 

{p,q)^ = j p{if^x)q{if^x)w{if^x)dx 

= / p{—i£^-\-x)q{—if-\-x)w{—if-\-x)dx 

Jr 

= 

Combined with Proposition this yields the asserted hermiticity property. □ 

We recall that the classical Hermite polynomials Hi (x) satisfy the orthogonality 
relation 


(19) 


f Hj{x)Hi{x)e ^ dx = 5ji2^l\\/TT, j,lGZ>o- 

Jr 


Combining this fact with the factorisation (14) of the intertwining operator 


it is now readily established by induction on the length n of A that the CEHPs 
H\.i{x) are orthogonal with respect to the Hermitian form (•,•) (cf. [T^l. 
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Theorem 1. The CEHPs H\ i satisfy the orthogonality relation 

71 

(20) =5jiV^2H\\{2{l-km), Z S Z>o \ {fci,..., fc„}. 


m —1 


Proof. The assertion clearly holds true for n = 0, with the empty product taken to 
be equal to one. Introducing the partition 


A — (A 2 ,..., A„), 


we have 


Ji^+S. 


Since Wm = the (formal) adjoint of Di is given by 

^ 1 ' 


d d 

Dl = - 

dx dx 


log 


W-2 


The factorisation 

thus entails that 


DlDi = .if; - 2fci - 1 


This completes the induction step, and the theorem is proved. 


□ 


Remark 1. Since (•,•) is Hermitian, each squared norm {p,p), p € Qa.Rj is real, 
but need not to be positive. In fact, if partition is not double, there is always 
a finite number of polynomials with negative squared norm, which can he easily 
identified using formula (20). For example, setting A = (1) in (2M, we see that 


(R(i)^i, R(i)_i) < 0 if and only if I = 0. Grinevich and Novikov 
similar fact in a finite-gap case. 


3j pointed out a 


We conclude this section by showing that the linear span U\ is dense in Qa,k in 
the sense that 

{p, Hx,i) = 0, V/ G Z>o \ {fci, ...,kn} p = 0. 

By Proposition we can formulate the result as follows. 

Theorem 2. The subspace Q\ is dense in Qa,r- 
Proof. Suppose that p G Qa,b is such that 

{p,q)=0, WqGQx. 

Introducing the polynomials 

q\,iiz) = Wx{z)Hi{z), I G Z>o, 
which clearly belong to the subspace Qx, we obtain 

0 ={p,q\,i)= [ p{z)Hi{z)e~''^dz, VZ G Z>o. 


Since the integrand is entire, we can take the limit ^ 0. Then expanding p in 

terms of the classical Hermite polynomials Hi, it follows immediately from (19) 
that p = 0. □ 
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Remark 2. If we assume that X is a double partition, then we recover orthogonality 
and completeness results from [T^] (see Propositions 5.7-5.8). Indeed, to recover 
the former it is enough to note that the weight function Q is guaranteed to he 
non-singular on the real line, so that we can take the limit ^ 0 in (18l; and the 

latter follows from the observation that we have Qa,k = C[z]. 


3. Exceptional Laurent orthogonal polynomials 

In this section we generalise our approach to the space of Laurent polynomials 
A = P,[z,z~^] using the trigonometric monodromy-free Schrodinger operators [4], 
which play an important role in the theory of Huygens’ principle [2]. 

More specifically, we consider the Laurent polynomials as defined in ([^. 

Due to the results of Theorem[^and Proposition]^ we call PK,a-u I & exceptional 
Laurent orthogonal polynomials (ELOPs). 


3.1. The general case. In this first subsection we allow any choice of complex 
parameters a = (oi,..., a„), € C \ {0}. 

We start from the elementary fact that the exponential functions 

ei(x) = exp{ilx), I G Z, 


have the eigenfunction property 

.ife/ = X G C/27rZ. 

Note that instead of usual unit circle M. 12711 we consider its complex version - 
cylinder P.12771. It is natural from the trivial mondromy point of view, see [1]. 

Sequences of Darboux transformations at the levels 0 < < A:„_i < ■ ■ ■ <k\ are 

now parameterised by n complex parameters B = (0i,..., 0„), 0^ G C. Specifically, 
introducing the functions 


(21) = 2cos(%a; + 0j), j = l,...,n, 
the resulting Schrodinger operator takes the form 

(22) ^ _ 2^ (log Wr(<)„„ ..., ^, J), 
where n = {ki ,..., A:„}. Furthermore, letting act by 

^ _ Wr((^,())fc,,...,(^fc„) 

Wr(0fc,,...,0fcj ’ 

the intertwining relation (|11[) holds true, and the functions 


(23) 




Wr(e;,())fci,...,())fc„) 


I G Z, 


have the eigenfunction property 



We note that at each level kj, j = 1,... ,n, the multiplicity is reduced from two to 
one. Indeed, by -([2^ and linearity of the Wronskian, we have the relation 


exp(i0j)</>K,fc^.(0j;x) + exp{-i9j)(j),^_k.{6j;x) = 0, j = 1,... ,n. 
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To establish the precise connection between the functions (j>K,,i and the ELOPs 
given by (|^, we change variable to 

2 = exp(ia;) 

and fix the values of the parameters a = (oi,..., a„) according to 




Ofe = exp(zdfc) 

£ C \ {0}, k = 1,..., n. 

Then, it 

is readily seen that 




(l>K,e-,i{9,x) 

= p.,aAz)n’.Az)~\ 

with P,^ 

ag{z) given 

by (j^ and 

‘^ki{ai;z) 

D$fci(ai; 2 ) 

^k2{a2;z) ••• $fc„(a„;z) 

D^k2ia2',z) ••• D^k„{a„;z) 

(24) 

W«.a{ 2 ) = 




-D”-^ $( 03 ( 02 ; 2 ) ••• 


where D = zd/dz and 

z) = az^ + a~^z~^. 

Furthermore, a direct computation reveals that PK,a\i is an eigenfunction of the 
operator 


= -D^ 


DW 

y y ^ y^K.a 


W,.a 


W. 


with eigenvalue P. 


Example 2. In the particular case k = {1} the corresponding Schrodinger operator 
is given by 

_ 2 

dx'^ cos2(a; + 0i) ’ 

When expressed in terms of the variable z and the parameter ai, the first few ex¬ 
ceptional Laurent polynomials P{i}^a-i defined by and the corresponding eigen¬ 
functions d){i} Q.; = ^ ^ given by 


P{l},a-,0 


aiz — aj^^z 


‘h{l},a;0 


— 1 —1 
aiz—a^ z 

aiz-i-a'^^z~^ ’ 


^{l},a;-l = 2ai, 

-P{l},a;l = —2ai 
P{l},a;-2 = O-l ^+3ai2: 
-P{l},a;2 = 01-2^ + ^Z, 


$ 


{l},a;-l 


2(21 


$ 


ai2+aj^ 2 ^ ' 

2a~^ 

ai2+a^^2“^ ' 
—1 —3,0 —1 

I a-, z +3aiz 


{l}.a;l 




ai2^+3aj^ ^2 
aiz+a^^ z ^ 


From these explicit formulae, it is manifest that both -P{i},a;±i 

linearly dependent and that each eigenfunction is singular at z = -F.i/a\. For general 

I G Tj, the latter fact can he easily seen from the definition o/-P[i},a;/- 
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We note that, upon setting 


J, 


the intertwining operator factorises according to (14|-(15). Just as in the 
Hermite case, it follows that each Schrodinger operator has trivial monodromy. 
Moreover, every monodromy-free trigonometric Schrodinger operator is of the form 
@, see g. 

Let Zk be the set of zeros Zi € C of the function WK^aiz) with multiplicities 
TOi € N and be the corresponding set consisting of Xj such that exp{ixj) = 
Zj, Zj G Zk (we drop the dependence on a in the notations for brevity in the rest 
of this section). 

Introduce the subspace 

= {P G A : $(a;) := (P/>VK)(exp(za:)) is quasi-invariant at all xj G X^}. 

It follows from trivial monodromy property that 

Wk := {P^,i : / G Z) C 

However, in contrast to Hermite case (see Proposition]^, the converse inclusion 
does not hold. Instead, we have the following result. 

Proposition 5. The codimension o/Uk in is n. 

Proof. From we deduce that 


Pm(^) =^'+'"'dety(Z,fci,... 


, kn) kj 
i=i 


l.d.. 


where 


\k\ 

i=l 

l.d. stands for terms of lower degree and V is the Vandermonde matrix 

^1 1 ••• 1 ^ 


V(Ol, ■ . ■ , 


Oi 


02 


\J^1 ^2 -m J 

Since det V{l,ki,..., kn) = 0 if and only if 1 = /ci, ..., fc„, it follows that the degree 
sequence 

It = {degP(z) :P GU^} 

stabilises at fci J- |k| J- I in the sense that I G for alH > fci J- |k| J- 1. Applying 
the same line of reasoning to the Laurent polynomials Pk,^-i{1/z), we find that the 
same statement holds true for 

C = {degP(z"^) : P GU^}. 

Among the ELOPs Pn,i with |Z| < fci J- |k| J-1, a maximal set of linearly independent 
Laurent polynomials is given by 

I G {ki,kn} U {0, ±1,..., ±{kn - 1)} U {±{kn -I- 1), ..., ±(fc„_i - 1)} 

U • • • U {±(^2 -I- 1),..., ±(fci — I)}. 


,m— 1 


a 


m — 1 


a, 


m —1 
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The cardinality of this index set equals 

Ti ‘2kn — It — kji — 1) T * * * T — ^2 — 1) “ 2/^1 — 71 -t- 1- 

Observing that 

2fci + 2\k\ + 1 — (2fci — n + 1) = 2\k\ + n, 
we conclude that the codimension of in A is 2\k\ + n. 

On the other hand, counting quasi-invariance conditions, we find that the codi¬ 
mension of in A equals 2\k\ and so the assertion follows. □ 


Remark 3. In contrast to the case of usual polynomials there are several definitions 
of the degree of a Laurent polynomial, but none of them are convenient for our 
purposes. Let us define the L-degree LdegP of a Laurent polynomial P = 
with Cp 0, Cg 0 as q if q > —p, and p if q < —p. If q = —p the L-degree is not 
well-defined since it could be both p and q. Under these assumptions 

LdogP^.i = \k\ -\-1, I G Z_|_ \ K, LdogP^.i = — |k| + I, —I G Z_|_ \ k, 

otherwise it is not well-defined. Note that the polynomials P^^ky and with 

undefined L-degrees are linearly dependent. 

Next, we consider a particular complex bilinear form on given by (§ with 
W = Wk, and establish corresponding Laurent orthogonality relations. A re¬ 
lated Fourier theory for more general algebro-geometric operators was studied by 
Grinevich and Novikov in m- 

Definition 2. Let p, G M>o be such that 

(25) Vz, ez,. 

Then, we define a complex bilinear form (•,•) on by setting 

(26) (P,Q) = ^/ p{z)Q{z)W-^-, P,Qg^,, 

271* Jc,, z 

where 

(27) Cp = {zGC-.\z\=pi}. 

Substituting z = exp(7a;) and following the line of reasoning used in the proof of 
Lemma 1^ we readily find that the product is well-defined in the sense that it does 
not depend on the choice of p,. More precisely, we have the following lemma. 


Lemma 2. For any P,Q G 
provided (251 is satisfied. 


the value of (P, Q) is independent of p G K>o 


We are now ready to state and prove the first of the main results in this section, 
which may be viewed as a natural analog of Theorem 


Theorem 3. The ELOPs P^y satisfy the Laurent orthogonality relation 

n 

=Sj+i,0 n 

m—l 

Proof. Just as in the proof of Theorem we note that the assertion holds true for 
n = 0, and proceed by induction on the length n of k. Letting k = (^ 2 , •.., fc„), we 
have 

1 

[PkP.Pk.i) = tt / {Di(j)k,j){x){Di(j)f,,i){x)dx. 

271 Jq 
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Making use of the factorisation 

( 28 ) DID^ kl 

with 



the (formal) adjoint of Di, we deduce 

which completes the induction step. □ 

Remark 4. Having started from an eigenvalue problem with doubly degenerate 
eigenvalues, we have that {Pk,,i, Pk,-i) — 0 for some of the ELOPs Pk,i- More 
specifieally, it is evident from the theorem that this is the case if and only if I = 
m = 1,... ,n. 

Expanding on the result of Proposition we proceed to establish the precise 
relationship between and We begin with a general definition. 

Let E be a vector space over C. Then V is called complex Euclidean space if it 
is equipped with a non-degenerate bilinear form B '.V ®V ^ E-. 

Definition 3 . Let W C V be a subspace of complex Euclidean space V. We say 
that V is a minimal complex Euclidean extension of W if 

dim (keri3|w) = codimylT. 

For any linear space W and bilinear form B with non-trivial kernel 

K ■= ker B, 

it is readily verified that there is a unique (up to isomorphisms) minimal complex 
Euclidean extension V D W. Letting K* denote the dual space of K, it can be 
realised as follows: 

V = K®K*® W/K, 
with the extension of B determined by 

(fci + ki+ Wi,k2 + k2 + W 2 ) i-t k2{ki) + ki{k2) + B{wi,W2), 

where ki,k 2 € K, fci, ^2 € K* and Wi,W 2 € W. Moreover, for each basis ki,... ,kn € 
K, there is a unique basis fci,..., G K* such that (kj, ki) = Sp. 

Example 3 . Suppose that B\yy =0, so that each vector w G W is isotropic. Then 
we have 

v^w®w*, 

with 

B{wi + Wi,W 2 + W 2 ) = W 2 {wi) + Wi{w 2 ), Wi,W 2 GW, Wi,W 2 GW*. 

As demonstrated by the following proposition, the inclusion Uk C provides 
a concrete example of a minimal complex Euclidean extension in the sense of Defi¬ 
nition |3] 

Proposition 6. is the minimal complex Euclidean extension ofU^. 
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Proof. From Theorem we infer that 

ker(-,-)lw„ = {PK,k, :j = l,...,n). 

(Note the linear relations OjPK.fcj + PK,-kj = 0.) Since ELOPs correspond¬ 
ing to different values of P, and hence different eigenvalues, are linearly indepen¬ 
dent, it follows that 

dim(ker(-,-)|wJ = n. 

Recalling Proposition we see that it remains only to verify that (•, •) is non¬ 
degenerate on Observing that 

Vjez, 

this follows, e.g., from the computation 

n 

= detV{l,ki,... ,kn)Y]_kj, 
i=i 

which is non-zero as long as I ^ dzkj , cf. the proof of Proposition □ 


3.2. The Hermitian case. In the case when all 9k are real or, equivalently, when 
parameters a = (oi,..., a„) satisfy 

(29) lofcl = 1, fc = 1,... ,n, 

we can introduce the Hermitian structure as follows. 

Note that in this case the weight function w{z) = 'Wk.{z)~'^ is invariant under 
the antilinear involution 


(30) 


P^(z) :=P(l/z), PgA, 


which will play much the same role as the Schwartz conjugate did in the Hermite 


case. In fact, observing that {DPy = —DP^ and that = ih/c, we can deduce 
from (24) that 

(31) 


wliz) = i-ir^--ppw,iz), = {fci,..., fc„}. 

In addition, the zero set Zf^ becomes invariant under the involution z 

Zi G Zf.^ y G Zf.^^ 


1/z, i.e. 


and, since z = Ifz whenever \z\ = 1, we have that 

(32) W,{z)Wl{z) = |z| = 1. 

Letting Z^ = {zi G Z^ : \z\ = 1} and X® = {xj : exp{ixj) = zj, Zj G Z^} C M, 
we introduce the following subspace of quasi-invariant Laurent polynomials: 

’^K,c = {P € A : $(ai) := (P/yVK)(exp(ix)) is quasi-invariant at all Xj G -A*} ■ 

From (HI), it is straightforward to infer that 

which allows us to define a natural sesquilinear product on ^k,c- 
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Definition 4. Assuming that (29) holds true, we introduce 

V = min \zi\, 

k.|>l 

and let /i G IR>o be such that 

(33) 1 < max(/i, l//i) < v. 

Then, we define a sesquilinear product on by setting 

(34) {P,Q)l = ^^£ P{z)Q\z){W^{z)wUz)y^^, P,Q e 

Again, the product does not depend on the specific choice of /i. 

Lemma 3. For any P,Qg ^k.,C: tbe value of (P, Q)l is independent of p G K>o 
provided (33) is satisfied. 

By adapting the proof of Proposition we can use the lemma to show that 
Definition yields a Hermitian product. 

Proposition 7. The sesquilinear product (•, •)l is Hermitian: 


{P,Q)l = {Q,P)l, yP,QG^.,c- 

Proof. Using the notation 

w{z) = l/W^,{z)Wl{z) 

and using a subscript to indicate the choice of p. in ( [M| ), we deduce the following 
equalities: 

Q)l,u = ^ / Pitie^nQiT-^eywitieydg, 


p27r 


= — / P{fie ^e^^)d(p 

27r/n 


= {Q, P)L,tj.-^, 

and so hermiticity follows from Lemma 

Moreover, the proof of Theorem is readily adapted to yield the following or¬ 
thogonality result. 


□ 


Theorem 4. Assuming that (29) holds true, the ELOPs P„_; satisfy the orthogo¬ 
nality relation 


(35) 


{Pn,j,Pn,i)L = Sji (fc^ - y, J, I e z. 


Proof. Taking z = exp(fa;) in the integral in (34) and observing that (cf. @) 
^(-x) = (-i)(”-™)(”-™+i)/^ir„(x), 
we establish the equalities 

1 _ 

{PK,j,Pii,i)L = -z- / {Di(j)f,y{x){Di(j)f„i){-x)dx 
27r Jg 


<.2-n- 


27r 


4>k,j{.x){DlDi(j)f,^i){-x)dx. 
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Appealing to the factorisation (28), we thus obtain the relation 
{Pk,],Pk,i)l = {kl - l‘^){Pk,j.Pk,l)L, 
and the assertion follows by induction on n. 


□ 


After replacing the bilinear form by a Hermitian sesquilinear form h, Definition 
l^as well as the succeeding discussion applies with minor changes also in the present 
situation. Specifically, we say that D is a minimal Hermitian extension of W if 

dim (kerli|vu) = codimyW. 

Then, we have the following analog of Theorem 

Theorem 5. The subspace which is the minimal Hermitian extension o/Uk, 
is dense in 

Proof. Suppose that P G ^k.,c is such that 

{P,Q)l = 0, 


Since the Laurent polynomials 

Q^,i = W,{z)Wliz)z\ IGZ, 
clearly are contained in we have that 

dz 


0 = {P,Q.,i)L = ^.f Piz)P — , VI gZ. 


Cu 


Taking the limit /i —?► 1 and using the property that 
1 

2717 

we conclude that P = 0. 


z^z^-=6k+i,o, fc,lGZ 
ld=i ^ 


□ 


Remark 5. It is known from the soliton theory that for every non-empty set k and 
any choice of real 9^ the corresponding potential always has singularities on the 
real line. This means that in the Laurent case we do not have non-trivial regular 
examples (unlike Hermite case with double partitions). 


4. Concluding remarks 

We have discussed two complex versions of the exceptional orthogonal polyno¬ 
mials, related to two classes of monodromy-free Schrodinger operators. We would 
like to emphasize two novelties compared to the original approach of Gomez-Ullate 

et al [TTl[T^ . 

First, in order to define the inner product in general we have to reduce the space 
of polynomials to the subspace of quasi-invariants, which has a finite codimension. 
The only exception is the Hermite case with double partitions considered in [12) . 

Second, in the Laurent case the space of quasi-invariants is not generated by the 
corresponding exceptional Laurent polynomials, so we need to consider the minimal 
complex Euclidean extension. 

In the rational case with sextic growth at inhnity there are some partial results 
[10) . which lead to finite sets of orthogonal polynomials of the same degree. It 
would be interesting to analyse this situation in a view of recent very interesting 
paper by Felder and Willwacher [3]. 
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It would be interesting also to see what happens with exceptional orthogonal 
polynomials in the multidimensional case. One can use the monodromy-free gener¬ 
alised Calogero-Moser operators, playing an important role in the theory of Huygens 
principle [5]. We plan to address this elsewhere soon. 
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